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1. Introduction 

There has been significant experimental progress owing to the remarkable success in B fac- 
tories. Recently there appeared measurements of the mass difference Am Bs from CDF |[T]] , Belle 
measurements of the pure leptonic decay B — > TV [§] , and the FCNC B — > p,coy . Also, sin(20! ) = 
«m(2j3)was measured with improved precisions. The semileptonic inclusive and exclusive decays 
b —>■ c,u were also measured with much higher accuracies. We can therefore overconstrain the 
CKM matrix elements with the present experimental data. This will be a good test for QCD calcu- 
lation , the standard model, and the physics beyond the standard model. 

The CP asymmetry Aqp{B — > J /yK) , the mass difference Am Bj j, the branching fraction of 
the pure leptonic decay SS(B —> T~v T ), and differential decay rates for various semileptonic B 
decays can be written as 

A CP (B -^J/yfK) « sin(20i) =sin(2j3) 

Am Bs = (known factors) m Bs f£ s B Bs \V ts V t b \, 
Amg, _ \Vt s \ 2 m Bs f B B Bs 
Am Brf \V td \ 2 m B Jl d B Bs ' 

^, s G%m B mz , m 2 . n7 
^(B^T-Vr) = F * T - 1 1 V| V ub Z X B , 

dT(B^D(*)lv) n c s |T7 l2 ( {w 2 -\) l l 2 F 2 {w) ForB^D* 
- (known factors) \V cb \ 2 ' 



dw ~ lVcm } (w 2 - \) 2 l 2 F 2 {w) For B - D 



d^T(B^X c lv) , l2 5r (function of Ai, AO 

.2., 1A 2 = ^ow 11 factors ) \Vcb\ 2 m 5 b [l + - 2 + • • ■ 

dt[dq z dm^ m L h 

dT(B -» 7T/v) G| 2 2 13/21,/ |2| f +/ 2ni2 

^2 = I ( v '^) "^1 l v «*l 1/ W )l . 



d 3 r(B^XJv) 2 5 (function of A 1; A 2 ) 
2 2 = (known factors) |Kfor«^[H 5 



The unitarity of the CKM matrix implies that |V /fe | = 1 + <9(A 4 ) and |V fi | = |V c6 |(l + <9(A 2 )), 
\Vub\ = |V c fc|A(p — iTj) and |V f <* | = |V c fc|A(l — p — ir\ + 0(X 2 )) using the Wolfenstein parameteri- 
zation. Thus, there are 11 independent experimental data for 3 unknown CKM parameters \V c b\, 
p and 77 in Wolfenstein parameterization. First, sin(20i) = sin(2j8), which is a function of p and 
tj can be determined purely from experiment. In order to determine the CKM parameters from 
other channels, we need to know the hadronic parameters: the decay constants f Bs , f Bd , the Bag 
parameters B Bd , B Bs , and the semileptonic form factors F, F*, and /+. HQET parameters nib, Ai, A2 
are also needed but they can be determined from experiment alone using the moments in inclusive 
semileptonic decays or rare decays. 

In fact, we already have strong constraints from Am B jAm Bd = IT .31^\\(stat) ±0.07 ps~ l 
, sin(20i) = sin(2/3) = 0.69 ± 0.03, and |y cfc | = [4.45 ± 0.045] x 10~ 3 with inclusive B -» X c lv 
decay. 

As can be seen from Fig. [I], the results are consistent with unitarity with 2 a level. However, it 
is also true that there is still large room for new physics. Since the error is dominated by theory ex- 
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Figure 1: Constraints on the unitarity triangle. The bands show bounds at 1 a level. 

cept for sin(20i) = sin(2/3), it is crucial to reduce the theoretical errors in the lattice determination 
of weak matrix elements for heavy flavor physics for more stringent tests of the standard model 
and the physics beyond. 

2. Heavy quark formalisms for heavy-light systems 

2.1 Lattice NRQCD 

Lattice NRQCD action is a discretized version of nonrelativistic effective action which is ap- 
plicable for the heavy quark mass whose spatial momenta are smaller than the mass. The expan- 
sion parameter is the velocity of the heavy quark for quarkonia and A/m for the heavy-light system 
where A is the typical momenta for all the light degrees of freedom. Since it is a nonrenormalizable 
theory, one cannot take the continuum limit. Also the action and operator can be matched to full 
QCD only by perturbation theory. In order to control the discretization errors the action is often 
highly improved at the tree level. The dominant source of errors are perturbative errors. 

2.2 Relativistic heavy quark (RHQ) formalism 

The Femilab action [Q] , AKT action [Q] , and the relativistic heavy quark(called RHQ) action 
[^] by RBC collaboration are the formalisms for heavy quark using improved Wilson fermion with 
suitably chosen improvement coefficients. The three formulations are essentially the same in the 
sense that they are Synamzik effective action applicable to quarks with small spatial momenta 
\ap\ <C 1 where the coefficients are mass dependent. These actions smoothly interpolate the static 
quark and light quark. Therefore one can in principle take the continuum limit without encountering 
the breakdown of the theory. However, since the discretization and perturbative error of the physical 
observable depend on am, how the B meson physical observable approach to the continuum limit 
is nontrivial. 

The discretization and perturbative errors are expected to be small by order estimation. Partial 
non perturbative (wavefunction) renormilzation using z'y° npert is useful. For higher accuracy both 
the discretization and perturbative error should be reduced. In order to reduce the perturbative error 
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either two-loop calculation which is possible only by automated procedure [Q] or nonperturbative 
renormalization. To reduce the discretization error further improvement by adding more terms is 
necessary. This is perused by FNAL, CP-PACS and RBC collaborations. 

Lin and Christ determined the coefficients of the RHQ action nonperturbatively in 

quenched QCD. 




They show that one can set ce = cb = cp by shifting ce and eg by field transformations 

W ^{l-a 2 [Y,f\[D\D^ E )w, (2-2) 
w ^{\-a 2 [f,f]{D\D^ B ) W , (2.3) 

so that only three parametersmo, £ , cp should be tuned. 

In order to determined the parameters nonperturbatively, they carry out the step scaling in three 
steps. In step 1, one starts with a very fine lattice in small volume on which am <C 1 is satisfied 





finer lattice 


coarser lattice 


L(fm) 


Size a }Lr 


sinze a-} arser (GeV) 


Step 1 


0.9 


24 3 x 48 5.4 GeV 


16 3 x 32 3.6 


Step 2 


1.3 


24 3 x 48 3.6 GeV 


16 3 x 32 2.4 


Step 3 


2.0 


24 3 x 48 2.4 GeV 


16 3 x 32 1.6 



Table 1: lattice setup for step scaling 

so that one can describe the heavy quark using Domain Wall fermion with controlled discretization 
error. One can then match the coefficients of the RHQ action on a coarser lattice for the same 
volume using one shell quantities: (1) the spin averaged IS state mass for heavy-heavy and heavy- 
light system, (2) hyperfine splitting for heavy-heavy and heavy-light system, (3) the spin-orbit 
average and splitting for heavy-heavy system, and (4) the dispersion relation. In step 2, 3 and so 
on, they can repeat similar procedure to match RHQ on a lattice to RHQ on an even coarser lattice. 
They demonstrate that one can actually determine the parameters with reasonable accuracy and 
obtain improvements in charmonium spectrum compared to those with perturbatively determined 
parameters. This method is quite similar to nonperturbative HQET by Alpha collaboration which 
will be explained later. However, at the moment the step scaling function is defined not in the 
continuum limit but a fixed lattice spacing assuming discretization error is under control. It will 
be important to have theoretical understanding about how the systematic errors in the matching 
procedure can be controlled in this method. 

2.3 Method with nonperturbative accuracy 

Rome II group [§J|,[0] proposed a method to compute B physics observables with nonperturba- 
tive accuracy based on finite size scaling. Consider a physical observable 0(Eh,Ei) which depends 
on two largely separated energy scales Ei and £/,(.£/ <C Eh). They assume that the finite size effects 
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has a mild dependence on high energy scale Ef,. Then Finite size effects can be obtained from the 
ratio Off of the observable in two different volume L and 2L. 

ff(E h E h ,2L) 

M£ " £ '" L) = ff(E h E h ,L) i2A) 
When Ei <C Eh the finite size correction can be expanded as 

o ff (E h E h ,L) = a e {E h L) + + + ■■■. (2.5) 

Eh E£ 

In the case of heavy-light meson almost at rest, the high energy quantity Ef, is the heavy quark 
mass and the assumption that one can expand the physical observable ratio in l/m is justified by 
HQET. Using the step scaling function Off one can obtain the physical observable in infinitely large 
volume as 

0{Ei,E h ,L„) = ^(Ei,E h ,L )Off(E,,E h ,L )o & (Ei,E h ,2Lo) ■ ■ ■ (2.6) 

When the volume is Lo small one can carry out lattice calculation with a cut off much larger than 
Eh with reasonable numerical cost so that one can compute the physical observable directly at 
energy scale Eh using the formalism of nonperturbatively 0(a)-improved Wilson fermion. But as 
the volume gets larger through step scaling at some point 2 k Lo becomes too large one cannot afford 
very small lattice spacing so that direct computation becomes hopeless. However, one can always 
find a lower energy scale E h < Eh where direct calculation is possible. In this case one can use 



Eq. 2.5 to extrapolate G&(Ei,E^ ,2*Lq) to Off(Ei,Eh,2 k Lo). Since each step can be extrapolated in 
the continuum with nonperturbatively 0(a)-improved Wilson fermion, the only systematic error in 
this procedure is the extrapolation in I /Eh. However, in order to take the continuum limit one has to 
know the parameter of constant physics so that one should know the master formula Xqcd scale and 
renormalization invariant quark mass as functions of the bare gauge coupling gf } and the bare quark 
mass mo. They find that in the case of the mass and the decay constant of the B meson mass one 
can practically control the extrapolation error at the level of few percent accuracy. The advantage 
is that this method is simple and promising. Probably, the Bag parameters, and form factors at zero 
recoil also fall into this category. Form factors for non zero recoil may be challenging. 

The Alpha collaboration proposes HQET with nonperturbative accuracy [ |l0| ] for high preci- 
sion computation in B physics. The action of HQET can be written with l/m expansion as follows 



L = L stat + ^ v >, (2.7) 

V 

(V)r(v) 



-(V) 

Utat = Wh [Do + Sm] ¥h , ^ (v) = I; ^ V) L\ V} (2.8) 
l| v ^ are the l/m v correction terms coj^ are their coefficients. 

= ¥h[-\o-B] ¥h , 4 1} = %[-^D 2 ] ¥h . (2.9) 

Since the static theory has a continuum limit and is a renormalizable theory, if we expand the 
l/m correction terms systematically to a fixed order n as operator insertions, one can renormalize 
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all the physical observable and take the continuum limit. In a very small volume where one can 
afford sufficiently fine lattice, the renormalization parameter can be determined nonperturbatively 
by carrying out QCD simulation for heavy quark with 0(a)-improved Wilson and imposing the 
matching condition for a set of physical observables {<J>jt(M,Lo)} as, 

tf QET (M,L ) = ^f D ( M ,L ) + O( ] ^),k = l,...,N„ (2.10) 

After matching QCD to HQET in small volume with lattice size Lo, one can then define the 
step scaling function Fk as 

^ QET (MaLo)=F k (^ QET (M,L )) + 0( ] ^ TT ),k = \,...,N n (2.11) 

By repeating this step, one can determine the matching conditions for (0- v ^ for large and coarse 
lattices where one wants to carry out lattice simulation. During each step one can take the contin- 
uum limit so that the only systematic error is the truncation error in 1 /M. To control the truncation 
error 1 /M <C Lo is required which restricts the smallest possible Lo as a function of M. This is 
in principle possible, but when one goes to higher order mixing with lower dimension operators 
through power divergences may give numerical difficulty, so that the calculation is technically more 
demanding. 



In this conference Guazzini et al. Ql lfl reported their proposal for further improvements. They 
combine the Rome II method and Alpha collaboration method. To be more precise, they basically 
follow the Rome II method, but hey also compute step scaling function a using nonperturbative 
HQET in the static limit. When they estimate the heavy quark mass dependence of the finite size 
correction, instead of extrapolating in 1/M, they make interpolation using the static result as an 
additional input. 



3. Heavy-light decay constants 

3.1 f Ds , f Bs in quenched QCD 

The determination of the heavy-light decay constants with nonperturbative accuracies is one 
of the most important progress. 

Since the charm quark is of order 1 GeV, the decay constant fo s in quenched approximation can 
be computed including nonperturbatively including the continuum limit with the present computer 
resources. Alpha collaboration JT2|]'s result for a -1 =2 — 4 GeV with 0(a)-improved Wilson 
fermion is 

f D ,= 252(9) MeV. (3.1) 

The Rome II group ^ computed the heavy-light decay constants in quenched QCD using 
0(a) -improved Wilson fermion by step scaling method. The observable is the nonpertubatively 
improved heavy-light axial vector current in SF boundary for vanishing boundary gauge field with 
periodic spatial boundary condition for fermions. They prepared three different size for L 3 x 2L 
volume with Lo = 0.4, L\ = 0.8,L2 = 1.6 fm for step scaling. The lattice spacings and RGI heavy 
quark masses are a = 0.012 - 0.033/m, m RGI = 1.6 - 7.0 GeV for L = L Q , a = 0.05 - 0.10/m, 
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m RGI = 2.0 - 3.5 for L = L\ and a = 0.10 - 0.20 fm, m RGI = 1.3 - 2.0 for L = L 2 . Defining the 
finite volume corrections factors with the ratio of the decay constants for two different volumes as 
o(Lq) = and cj(Li) = j^2L^) ' tne decay constant in the infinite volume can be obtained as 

= /u,(Io)<T(Io)cy(Li). (3.2) 

The result is 

f Bs {Lo) = 475(2)MeV, f Ds (L ) = 644(3)MeV (3.3) 

GbXLx) =0.417(3), Ob^Lj) =0.414(3) (3.4) 

a Bs (U) =0.97(3), o Ds (U) =0.90(2). (3.5) 

As it turned out, the heavy quark mass dependence of the step scaling function are indeed small, 
which justified the extrapolation. Combining these results 

f B , = 192(6)(4)MeV, f Ds = 240(5)(5)MeV. (3.6) 

Alpha collaboration compute static heavy-light decay constant with lattice HQET which is 
matched to QCD with nonperturbative accuracy by Schrodinger functional method. They computed 
the renormalization group invariant matrix element which can be related to the decay constant 



by a matching factor Cps [14] as &rqj = fps^ m Ps/Cps and obtain 

<D^ = 1.74(13) (3.7) 




Figure 2: Interpolation of static and relativistic results of heavy-light decay constant to obtain fs s . Figure 
taken from [ p3[ ] . 



Alpha collaboration [ |I5| ] also computed the decay constants for the charm quark mass regime, 
i.e. rag = 1 .7 — 2.6 GeV, at four lattice spacings in the range a = 0.05 — 0. 1 fm using 0(a) -improved 
Wilson fermion for both the heavy and the light quarks. They then interpolated the decay constants 
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in the static limit and those for finite quark mass to obtain f Bs . They found that both linear and 
quadratic interpolations lead to 

f B , =206(10)MeV. (3.8) 
using ro = 0.5 fm for the scale input, as shown in Fig.2. 



Guazzini et al. [11] reported a quenched study of f Bs with nonperturbative accuracy. Their 
approach uses the combination of the two methods. They compute the heavy-light decay constants 
in finite volumes both for the relativistic and static heavy quarks to the step scaling and "interpolate" 
the finite volume correction f Bs to using both the relativistic and the static. They computed the 
2 point functions for static and relativistic heavy-light axial current with Schrodinger boundary 
conditions, where the boundary gauge fields C = C = and periodic boundary condition in the 
spatial directions 6 = for the light quark. The data for relativistic heavy-light current is obtained 
by the reanalysis of those by Rome II collaboration []|]. They chose fhi^/mhi for the physical 
observable rather than //,/. Thus finite size corrections 0\, 02 are defined by the ratio of fhi-Jmu 



for different volumes as 



= fhi{2L {) )^/m^(2L Q ) _ fhl {L^^m(u) n „ 



the infinite volume can be obtained as 



fB s { L °°)\l m B s {L°°) = fis s (L ) y m Bs (L )o\Oi. (3.10) 

As shown in Figs. 3, the heavy quark mass dependences of the finite size corrections have much 
better control with the help of static results. Their preliminary quenched result is 

f Bs = 1 86 ± 6 MeV from Static + Rome II (3.11) 
f B , = 195 ± 1 1 MeV from only Rome II (3.12) 

which are consistent with previous results by Rome II and by Alpha collaborations. 

There are also calculations of heavy-light decay constants with Ginsparg-Wilson fermions. 
The RBC collaboration [ jl6[ ] has carried out a quenched study of D meson using domain wall 
fermion and DBW2 gauge action. The quark mass ranges from m q = \m s ~ \m s and the lattice 
spacing is a ~ 3 GeV. Using the nonperturbative renormalization factor for the light-light axial 
vector current Jl7| ] and giving the mass correction as 

7 hl rjlljionpert ^q,DWF{ am heavy) ,~ n , 

z a — z a 7 r> (J.ijj 

£q,DWF{amiight) 

their result is 

f D , = 254(4) (12) MeV, (3.14) 



where the errors are statistical, and systematic errors. Chiu et al. Ql 8Q, [[19|] also computed fjj in 
quenched QCD using the optimal domain-wall fermion on a lattice with a" 1 = 2.2(GeV) for 30 
quark masses am q = 0.03 — 0.80 using f n as scale input to find 

f Ds = 266(10) (18) MeV, (3.15) 

where the errors are statistical, and systematic errors. 
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o\ Decay constant: Static + Step scaling 



Plot of a t for the meson decay constant, static + Romell 

Without static If: 1.194(13} 
With static If: 1 .2050(88) 




With static qf 

With static It 
Without static 



0.04 0.06 0.08 

1/(L *M PS (L.)) 



-1.194(13) without the static point 
-1.2050(88) static point included 



0.12 0.14 

1 m PS^l'' 

0(L 00 ) =O(L a )a(2L )a(L xl ) 



(72 Decay constant: Static + Step scaling 



Plot of o for the meson decay constant, static + Romell 



Without static If: 1.002(54) 



With static If: 0.935(29) 




— With static qf 
With static If 
Without static 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 

1/(2L 1 'M ps (2L 1 )) 



erg = 1.002 ± 0.054 without static 



O(Loo) =O(L )a(2L )a(L oo ) 



0.935 ± 0.029 static included 



Figure 3: 1/M interpolation of the finite size corrections 0"i (top) and a% (bottom) for fg s Jmb s - Figures from 
Guazzini's talk. 



3.2 fo s , /b, in unquenched QCD 

FNAL/MILC collaboration [ ^0| ] reported preliminary results of fs s for iif = 2 + 1 flavor QCD 
with MILC configuration. They use fermilab formalism for the heavy quark and improved stag- 
gered for the light quark. The lattice spacings are a = 0.090.12,0.15 fm. The renormalization 
factor Za is taken to be 

= PfV^Zf, (3.16) 
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where Zy's are computed nonperturbatively and the remaining part is computed by one-loop 
perturbation theory. Their preliminary result is 



f B , =253(7)(41) MeV, f B Jf Ds = 0.99(2) (6), 
where the errors are statistical error and systematic errors. 



(3.17) 



*(/*)/*„, 




SF scheme, N,=2 

1— loop 7, 2-loop (i 

2- loop 7, 3-loop i 



100 1000 

m/a 



Figure 4: Step scaling of heavy-light axial current. Figure provided by J. Heitger. 



Alpha collaboration [21] computed the renormalization factor for the static heavy-light axial 
vector current Z% at for rif = 2 unquenched QCD. Their preliminary result is shown in Fig. 4. Their 
preliminary result is 



$>{L max )/<S> RGI = \.U{l), 



(3.18) 



where L max is the physical lattice size in which one wants to carry out the matrix element calcu- 
lation. Using this result, once the large volume rif = 2 unquenched calculation j8 = 5.3 for the 
regularization dependent renomalization factor Zf at (L ma x,go) and the lattice bare matrix element 
fg tat Jm.B s ) lat {Lmax , go ) is done, one can obtain the static heavy-light decay constant as 



<$>RGI 

{Lmax 



■Z^ (L max , go) (Ib"' V™^)'"' (Lmax, go), 



(3.19) 



where Cps is perturbatively calculable conversion factor. The large volume rif = 2 simulation is 
now in progress for j8 =5.3. 

3.3 Discussion on fs s , /d s results 



Fig. 5 show the summary of decay constants /b,,/d s in quenched, rif = 2 unquenched, and 
nf ■ = 2 + 1 unquenched lattice QCD. It should be noted that the quenched results are getting very 
precise owing to the recent developments with finite volume technique which allows us to compute 
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N=2 + l 

HPiQCD^OpS), 
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f B [MeV] 



Figure 5: Decay constants fojleft), /^(right) 



0.7 




Alpha 
■ Romell 
♦ Rome II + static 
FNAL 

Collins et al.(NRQCD) 
JLQCD(NRQCD) 



0.3 - 



0.2 



0.0 0.1 0.2 0.3 0.4 0.5 

l/m H (GeV~') 



Figure 6: Comparison of 1/m dependence of f\[m 



nonperturbatively renormalized heavy-light decay constants in the continuum limit as discussed in 
the previous subsection. We take the average of the results from Rome II and Alpha collaboration 
as the best result in quenched approximation, 



n f =0 
JR. 



194(6) MeV, f%~° = 245(6) MeV, ( 

\JD S 



n f =Q 



: 0.80(6), 



In the unquenched case, the decay constants have larger errors from perturbative matching. I would 
quote the average of HPQCD/MILC and FNAL/MILC results for f B and FNAL/MILC results for 
/o as the best value. However, since the best result come from the same configuration, it would be 
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worthwhile to study the heavy quark mass dependence and n f dependence based on the collection 
of results from various collaborations. 



n f Group 







MILC [|2J] 
FNAL [g] 
Lellouch-Lin [25] 
Rome II [Pi 



Alpha Q15p 
RomeII+ Alpha [|l|] 



Wilson 

fermilab 

clover 

clover 

clover 

clover 



Wilson 

clover 

clover 

clover 

clover 

clover 



heavy light 5^ 



0.89(4) 
0.88(3) 
0.82(5) 
0.80(4) 
0.81(6) 
0.76(3) 



2 MILC [||] 


Wilson 


Wilson 


0.92(7) 


2+1 FNAL/MILC [g§] 


fermilab 


Imp Stag 


0.99(2)(6) 



Table 2: The decay constant ratio fgj '/d s 



Table 2 shows the ratio of /bJ/d,- Recent quenched calculations show smaller values of jjf-. 
Fig. 6 shows the comparison of l/M dependence of fs s \J fm Bs m quenched QCD near the static 
limit by Alpha, Rome II, FNAL, Collin's et al. and JLQCD. It can be seen that the 1 /M slope is 
consistently small independent of the action or collaborations. Parameterizing 



fB s J fm Bs — {fB s J fm Bs ) itat ( 1 



_£L 



(3.20) 



both the Alpha collaboration and Collins et al give the slope ofci ~ 0.5 — 0.6 GeV. 

MILC results for rif = 2 suggest that sea quark effects may increase j& but not significantly 
due to the error. On the other hand, FNAL/MILC preliminary rif ■ = 2 + 1 result presented in this 
conference suggests a significant increase in the jjf- . However, one should bear in mind that the the 
systematic error is slightly different for B and D in fermilab formalism so that some consistency 
check is desired. 



Group 



heavy 



a 1 input 



■T 



•"/ 



■fo. 



n f =0 



73m 



jif=a 



JLQCD [27], [28] 



11 1 r 



CP-PACS [£ 
HPQCD [|3C 



NRQCD 
NRQCD a 
NRQCD r 



1.13(5) 
1.10(5) 



1.15 



CP-PACS [0] 
MILC pi, [Ee 



fermilab 
Wilson 



nip 



1.14(5) 
1.09(5) 



1.07(5) 
1.08(5) 



Table 3: n/ dependence of fs s , fo s - 



Table 3 is the collection of the rif dependence of the heavy Tight decay constants fs s , /d s using 
the same gauge and fermion action by the same group. It is seen that if the scale is set by the low 
energy inputs, turning on the sea quark effects from rif = to rif = n + 2 to increases /# t by 10-15% 
, while the increase is not significant for /b t . It is quite natural to expect the size of the sea quark 
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effects for fy s should be something between that for fg s and fa. And with the low energy inputs 
fx receives almost no sea quark effects by definition, the sea quark effects for fy s should be larger 
than for fo s , which explains the above observations. 

From Tables [O], [3.3[ , we also estimate the ratio of decay constants as 

f f=1 f "f= 2 

= 1.12(5), ^o = 1.08(5), (3.21) 

Jb, fo » 

r"f= 2 + l a/ =2+l 

■^T = L15 ( 5 )' tjh ^r = 1-10(5). (3.22) 

This can give an educated guess for tit = 2 decay constants. However, there are several uncertain- 
ties in this argument. First, the up/down sea quark mass for the unquenched configuration other 
than MILC may not small enough to fully reproduce the sea quark effects. Also when one uses low 
energy inputs the scale suffer from the chiral extrapolation uncertainty. Although Sommer scale 
r is relatively stable, but the phenomenological value of r = 0.5 fm also suffer from uncertainty 
which is typically 10%. Our educated guess for rif = 2 results are 

f n B '= 2 = 217(12)(22) MeV, f D [ =2 = 265(14) (27) MeV, (3.23) 
f B [ =2+l = 223(17)(22) MeV, f^ 2 = 270(18)(27) MeV, (3.24) 

where the second error is added to take account the scale uncertainties of order 10%. On the other 
hand the average based on the actual data of decay constant in nj ■ = 2 + 1 QCD by FNAL/MILC 
and FNAL/MILC collaborations are 

f n B [ =2+l = 260(30) MeV, /g =2+1 = 249(16) MeV, (3.25) 

which is marginally consistent with our estimate within errors. Combining my educated guess and 
HPQCD/MILC, FNAL/MILC results my 'world average' would be 

f n B [ =2+l = 240(30) MeV, /^ =2+1 = 260(20) MeV. (3.26) 

3.4 chiral extrapolation 

In order to obtain fs d and fo d one has to take the chiral extrapolation. This offers another im- 
portant issue for precise determination of the decay constant in addition to the problems discussed 
for fg s and fo d . The correct answer can only be obtained with unquenched calculation. The chiral 
perturbation theory tells us that the chiral logarithmic corrections to the SU(3) breaking ratio of the 
decay constants is fl33"l] 



fB sy /mB s _ , , 1 + 3g 2 
fB d JmE 



FNAL/MILC collaboration pOp reported preliminary results from nf = 2+1 heavy-light de- 
cay constants in the previous subsection. With the staggered quark the pseudoscalar mesons for 
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each flavor quantum number (7) has 16 tastes labeled by E, = P,A, T,V,I. Their masses are splitted 



as 



2 



j = {m a + m h )ii + a l A^ %=P,A,T,V,I (3.28) 



The staggered chiral perturbation theory suggests that the quark mass dependence of the heavy- 
light decay constant 0> H? = f HgV /mJ^ is 

®H q = <s>H [i + + anal y tic terms ] . ( 3 - 29 ) 

where the explicit form of 8/h , which is the analog of the chiral log in the continuum chiral 
perturbation theory, can be obtain from staggered chiral perturbation theory [f35|]. Due to the taste 
symmetry breaking of 0{a 2 ) terms they have many parameters for a 2 effects which have to be fitted 
from the lattice spacing dependence of the lattice data. Some parameters can be obtained from the 
pion system but other parameters have to be fitted from the data of the heavy-light decay constants 
themselves. Their preliminary results are 

(f D JfD d ) nf = 2+l = 1-21(1)(4). (f B JfB d P= w = L27(2)(6), (3.30) 

where the errors are statistical and systematic errors. 

Gadiyak and Loktik [ 32 ] made a rif = 2 unquenched study of SU(3) breaking effect using 
domain wall fermion and DBW2 gauge action. The quark mass ranges from m K = 490,610,700 
MeV and the lattice spacing is a ~ 1.69(5) GeV. They found that 

(fB s /f Bd r= 2 = 1.29(4)(4)(2). (3.31) 



Group 



heavy light 



"/ fB s /fB d 



visible chiral lot 



CP-PACSS 
CP-PACS |fl) 
MILC \M\ 



JLQCD [28] 



Gadiyak and Loktik [E 



NQCD clover 2 

fermilab clover 2 

fermilab Wilson 2 

NRQCD clover 2 

static DW 2 



1.18(2)(2) 
1.20(3)(3)(12) 



NO 
NO 



1.16(1)(2)(2)(^) NO 
NO 
NO 



1.13(3)0 
1.29(4X6) 



HPQCD/MILC 
FNAL/MILCC Q20I] 



NRQCD Imp Stag 2+1 1.20(3)(1) 
fermilab Imp Stag 2+1 1.27(2X6) 



YES 
YES 



Table 4: SU(3) breaking ratio f B Jf Bd 



Tables 4,5 show the collections of the unquenched results of fuj fD d and fyJfB d - Except 
for FNAL/MILC and HPQCD/MILC, they do not observe the chiral log. This is natural because 
other results use much heavier light quarks. Fig.7 show the comparison of the light quark mass 
dependence of fg s <JmB s / fB d \/™B~d fr° m JLQCD and HPQCD. They show consistent behavior for 
larger light quark mass. It seems that the JLQCD result may be missing the possible onset of chiral 
log which is found by HPQCD data. However, the results with MILC configuration are obtained 
through the staggered chiral perturbation theory, which requires quite complicated analysis with 
many parameters. Independent calculations with other formalisms are needed. 
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Group heavy light rif /dJ/dj visible chiral log 


CP-PACSS [ 
MILC [|2§] 


31 


] fermilab clover 2 1.18(4)(3)(^) NO 
fermilab Wilson 2 1-14(1)(+|) (2)(+g) NO 


FNAL/MILC [ 


20 


] fermilab Imp Stag 2+1 1.21(1)(4) YES 



Table 5: SU(3) breaking ratio f D Jf Dd 



• JLQCD (2003) 
■ HPQCD (2004) 




0.0 



0.5 1.0 

m /m 

q s 



1.5 



Figure 7: light quark mass dependence of &B s /&B q , where 4> = ^/msfs ■ 

4. Bag parameters 

The bag parameters that parameterizes the B® — B® mixing amplitude are defined by 





-YsW\B q ) 


= \m\f B B Bq 


( where q = d,s), 


(4.1) 


<B^(1-75V^(1 


-Ysy\B° s ) 




( where / ? = (m " fo + ^ ) ), 


(4-2) 


(B° s \b\\-Y 5 )qJbJ(l 




_ 1 m 2 f 2 Bs 
- 3 m B q JB q W ■ 




(4.3) 



HPQCD J3q ] computed the bag parameters for B s mixing calculation with improved rif = 2 + 1 
dynamical staggered quark. The simulation was carried out using NRQCD action for heavy quark 
and AsqTad action for light quark on a 20 3 x 64 lattice with a -1 ~ 1.6 GeV with the valence light 
quark mass &Xm s and the ud sea quark mass at 0.25/n s , 0.5m v . They computed the matrix elements 
for three types of AB = 2 four fermion operators which correspond to f^Bs,, /| s ff an d /| §f ■ 

Defining AB = 2 four-fermion operators as 



OL = [b l q'}v- A [b J q J }v-A, 
OS = [&y] S -H&V]s-P. 



(4.4) 
(4.5) 
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<23 = [LVh-pftfis-P, (4-6) 
OLjl = ^[W'r^v-A^^lv-A + [Vj]v-A\ybi?q>]v-A, (4-7) 

OS/1 = ^_ [VP'fA-p^is-p + PV]5-p[VP7^']5-p, (4-8) 

23 jl = [V^'7^] 5 _p[P'^] 5 _ P + [frVls-plW^j-P, (4-9) 

where z, j are color indices. The former three are operators in the leading order in 1 /M and the latter 
three operators are 0(1 /M) operators. The lattice operator is matched to that in the continuum 
using one-loop perturbation theory as 

1 (OL) Ms 



2m Bs 
1 

2m Bs 
1 

2m Bs 



(OS) MS 
(Q3f S 



+ [1 + ap LL }(OL) eff + apLs(OS) eff 


(4.10) 


+[(OL) eff - a{ffi(OL) eff + Cw(OS) ef f)], 


(4.11) 


+ [1 + ap ss ] {OS) eff + aps L (OL) eff 


(4.12) 


+{{OSj\) eff - a(ffi(OL) eff + Cfi(OS) eff )], 


(4.13) 


+ [1 + ap^{Q5) e ff + ap3 L {OL) eff 


(4-14) 


+[{&fl)eff ~ a(tf L (OL) eff + Ci 3 o 3 (G3) e// )]. 


(4.15) 



It should be noted that in this work dimension 7 operators are included for the first time in NRQCD. 



Previous work by Hiroshima group p8| ] and JLQCD [|39J, [28|] include only dimension 6 operators. 

They find that the sea quark mass is only a few percent and quote the result for m sea = 0.25ra. v 
as their best value. 

Jb s ^B\ s = 0.281 (20)GeV, / BiV /fi^ = 0.227(16)GeV, (4.16) 
/ Bs W^ = 0.295(21)GeV, fay/k = 0.305(21)GeV (4.17) 

The key points of HPQCD's result is that the direct calculation of f B B Bs gives better accuracy 
than computing f Bs and B Bs separately. The bag parameter has a smaller central value than that 
of JLQCD (rif = 2) after including 1/M correction (dime=7 operator), which is not included in 
JLQCD's calculation. On the other hand, f Bs has a larger central value than JLQCD so that f B B Bs 
is consist ent T able 6. 

/fl t \j B Bs is related to the mass difference in B s — B s mixing as 

Am B , = -S- m B JlB B m 2 w S (mj/ml)\V ts V tb \, (4.18) 

n b 

where 77 is perturbatively calculable factor and So(mf /m^) is the Inami-Lim function. CDF [[]]] 
recently measured the mass difference as 

Am Bs = 18.3(^)ps _1 (4.19) 
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Combining this with recent \V c b\ value and CKM unitarity relation, the above equation requires 

^7^ = 0.245(20) (GeV). 

Using heavy quark expansion the width difference of B s — B s mixing can be obtained at NLO 

as 



AT\ \6k 2 B{B s Xev) J%m B , 



vg 2 



r J b s 8{m 2 c /m 2 h )f]QCD m\ 

x (g(z) \b b , (m h ) + G s (z) 5 - + ^/l-W/^5i/ m ) , (4-20) 

where g(z) = 1 — 8z + 8z 3 — z 4 — 12z 2 lnz and tjqcd is the short distance QCD correction. G(z) 
and Gs{z) are NLO QCD corrections which appear in OPE. 8i/ m is the NLO contribution in I /mt, 
expansion. Using HPQCD results they predict 



(jr\ =0.16(3)(2), (4.21) 



where the second errors comes from the uncertainty in the correction term wl —4mc 2 /m^di/ m . 



n f 


group 


heavy 


B B s ( m b) 


B Bs {m h ) 


BB s (m h ) 





Becirevic. et al. [|37|] 


HQET 


0.87(5) 


0.84(4) 


0.91(8) 





JLQCD @ 


NRQCD 


0.84(5) 


0.85(5) 




2 


JLQCD [||] 


NRQCD 


0.85(6) 






2+1 


HPQCD [||] 


NRQCD 


0.76(11) 


0.84(12) 


0.90(13) 



Table 6: The bag parameters and Bg (m&) 



n f group 



heavy B Bd (m b ) 



UKQCD 

Becirevic. et al. [37] 

JLQCD @ 

2 JLQCD @ 



HQET 0.87(5) 

HQET 0.87(6) 

NRQCD 0.84(6) 

NRQCD 0.84(6) 



Table 7: The bag parameters and B^(mi,) 



Table 7 gives the summary of B Bd from various collaborations. It should be noted that HPQCD's 
result with 1 jm correction in the operator gives lower values. Further understanding of 1/m depen- 
dence is required. On the other hand, the light quark mass dependence seems small from the data. 



In fact, chiral perturbation theory [33] suggests that the light quark mass dependence is 



B B , _ 1 



— ml\ogml + - 



(4.22) 



Since g ~ 0.6, The coefficient of the chiral log is very small, which agrees with the lattice results. 
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5. B — > nlv form factors 

The matrix element (■K{k K )\qy il b\B(pB)) for the heavy-to-light semi-leptonic decay B —* nlv 
is often parameterized as 



{n{k 7l )\qfb\B{p B ))=f + {q 2 ) 



(pB+k.V-^-^q" 



+ f{q2) n^ Lq ^ (51) 
1 



with ps and ^ the momenta and q = p B — k % . The differential decay rate of the semileptonic 
B° — > K~l + Vi decay is 



1 G 2 



^^I 3 |/ + (^ 2 )| 2 - (5-2) 



|Kfc| 2 ^ 2 247T 3 ' 

from which one can extract the CKM element \V u b\. 

HPQCD collaboration p0| ] has made a new study of B — > 7i/v form factors using 2+1 flavor 
MILC configuration with t? -1 = 1.6,2.3 GeV. They used NQCD action for the heavy quark and 
improved staggered fermion for the light quark. The light quark mass ranges m q /m s = 0.125 — 0.5 
on the coarse lattice and m q /m s = 0.2 — 0.4 on the fine lattice. The heavy-light vector current 
is renormalized with 1-loop matching through 0{a/M) and 0{aa). The chiral extrapolation is 
carried out using staggered chiral perturbation theory. In order to make the analysis convenient 
they parameterize the matrix element as 

(K(k„)\V^\B{p B )) = s/2^[v^f }l +k^f ± ], (5.3) 

with 



El, ^ = *£-(Vv)v". (5.4) 



In order to interpolate in q 2 , they used several different pole model fit functions. The first one is 
BK parameterization with three parameters with q 2 = q 2 /mg* , 



/ + fa 2 >= ( ,-^-^ r /V) = ^. (5.5) 



The second one is BZ parameterization with four parameters 



/ + (<? 2 ) = 7T%+ n Sr -p* , (5-6) 
The third one is RH parameterization with four parameters 

M-fm^LfL. ( 5.7) 

(i-r)(i-r/r) 

For all three cases the parameterization is the same for f°, First, the momentum dependent form 
factor data is interpolated to fixed E^s using these parameterization. Second, the chiral limit is 



taken for each fixed E n using staggered chiral perturbation theory [35]. It turns The results with 
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n f 


Group 






heavy 


r(r>i6G^) i 








UKQCD 
APE [15 
FNAL [4 
JLQCD [ 


i 

6] 
47 


i 

i 


clover 
clover 
fermilab 
NRQCD 


2.30(+^)(51) 
1.80(^)(47) 
1.91(1^X31) 
1.71(66)(46) 


2+1 
2+1 


HPQCD/MILC @] 
FNAL/MILC [0] 


NRQCD 
fermilab 


1.46(23)(27) 
1.83(50) 



Table 8: Values for partial branching fraction [c ' , |i e — ' ps 1 for various lattice calculations. 



different parameterizations are very well consistent with each other. Choosing BZ parameterization 
for the best result, they obtain 



1 rii ax dT 



\V u b\ 2 Jl6GeV 2 dq 2 

Using the experimental data from HFAG @] Br il 2 > 16GW 2 ) = 0.40(4) (4) x 10~ 4 , it leads to 

\Vubl =4.22(30)(51) x KT 3 , (5.8) 



which should be compared with FNAL/MILC results [|41|], 

\V ub \ = 3.76(25)(65) x 10~ 3 (5.9) 

(See Fig. 8). Table || shows the partial branching fraction for q 2 > \6GeV 2 for various lattice 
calculations. So far within large errors, all the results are consistent. The average of nf = 2 + 1 
results seems somewhat smaller than that of nf = but not significantly. 

6. mi, 

Alpha collaboration made a quenched study of l/M corrections to HQET, which is an update 
of last years work. Matching of QCD and HQET at small volume, step scaling in HQET towards 
larger volume and computation of ms s in large volume and finally convert to nib. Last year they 
quoted that value 

mf = mr+m^\ (6.1) 
mf at = 4.350(64)GeV, = 0.049(29)GeV (6.2) 

Guazzini et al. [|ll|] computed the bottom quark mass using similar finite size scaling as fs- 
Their preliminary results are 

mf^ = 6.96(11) GeV (only Rome II), (6.3) 
ml G1 = 6.89(11) GeV (Static + Rome II). (6.4) 
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Figure 8: B — > nlv form factors in rif = 2 + 1 QCD. 



Kronfeld and Simone p9| ] made a quenched study of HQET parameter A, Ai, and m^. The 
idea is that HQET relation the heavy-light meson mass can be expressed as 



M(m) = m + A + 



2m 



dj^ + 0{\/m 2 
2m 



(6.5) 



Fitting the mass dependence of the heavy-light meson from lattice calculation one can extract A, 
X\ , and mi, in lattice scheme. Using perturbation theory one can then convert HQET parameters to 
another short distance scheme free from renormalon ambiguities. Application of this method by to 

nj 



2 + 1 unquenched QCD by Fermilab collaboration is reported in this conference [|50|]. 



n f Renormalization Group 



mb(mb) (GeV) 



ms and HQET 







NNLO 
Non pert. 



Martinelli, Sacrahjda [0] 4.38(5)(10) 
Delia Morte et al. |E2] 4.350(64) 



2 NNNLO 
NNLO 



Renzo et al. [53] 



McNeile et al. [p 



4.21(3)(5)(4) 
4.25(2)(10) 



T and NRQCD 



2+1 



NLO 
NLO 



Gray et al. 
Nobes, Trottier 



4.4(3) 
4.7(4) 



Table 9: m b (m h ) 



Recently HQET parameters are extracted by the global fit of the various moments for inclusive 
B decays such as {Ef), (mj?) in B -> XJv or (Eft, in B -> X s y . The result [|^] is 

mp = 4.20(2) (5) GeV, 
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where the first and the second errors are the experimental and theoretical errors. These deter- 
mination is used to improve the accuracy of \V c b\ and \V u b\ determinations from inclusive semilep- 
tonic decays. A better determination of HQET parameters would provide further improvement in 
\V c b\ and \V u b\ determinations, which will be possible near future. Summary of recent results are 
given in Table 9. 

7. New methods 



7.1 Dispersive bounds for form factors 

The momentum range of B —> nlv form factors computed from Lattice QCD is limited by the 
small recoil or large q 2 region. This leads to a big disadvantage because most of the experimental 
data lies in large recoil region. While one can extrapolate in q 2 with a fit ansatz, this will always 
introduce some model dependence. Dispersive bounds is one possible way to constrain the q 2 
dependence in model independent fashion using unitarity. 

IL \L 



Figure 9: A map from t plane to z plane 



Consider the imaginary part of the vacuum polarization amplitude for the current V(x) 
uy^b{x) and a map as in Fig. 9 



U^iq) = iJd 4 xe^ x {0\T{V^x)V vJf (0)}\0) 
z(t,t ) = ^^-^I=p< ( with t ± = (m B ± m7[ ) 2 ), (7.2) 



Then, from dispersion relations one obtains 



, is 1 d 2 r ~ n 1 f°° , thalliU 



2d{q 2 ) 2 ^ iJ nk {t + Q 2 f 

with Q 2 = —q 2 and Y\ an isospin factor, while % 's can be computed using OPE and perturbative 
QCD. Unitarity tells us that this is equal to the sum over all the hadronic states, and dropping all 
the excited states and leaving only the B % state gives an exact bound. 



3 



487T t 

^ K '- t ^-'- )1 'W <MW. (7-4) 
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shows that an upper bound on the norm can be established by choosing [recall that \z\ = 1 along 



the integration contour in ( |7.3| )] 

MMo) = 
00 (Mo) = 



t+-t f z (t,0)\fz(t,-Q 2 )\ 3/2 fz(t,t )\ --fz[t.t- 



4SK(t+-t y/ 4 



Q 2 -t 



t Q -t 



-1/2 



-3/4 



t?m_ V^W (z(t,0)\ ( z(t,-Q 2 ) \ /z(Mo) V 1/2 
16tt (^-^o) 1 / 4 V -f )\-Q 2 -t )\to-t 



-1/4 



.(7.5) 



Combining Eqs. [Lp , [7.4 |775| and making change of variables in the integration from f to z. 
We obtain 



<0/o|0/o) < Xo, (PM+\PM+) < X+, 



(7.6) 



where J is a quantity which can be obtained using OPE and perturbative QCD. The inner product 
(g\h) for arbitrary functions g(z) and/i(z) is defined by the integral along the unit circle in z plane 
as 

ch , , (7.7) 



^*(z) = z{t,m* B ) is multiplied to / + in order to remove ZT pole inside the unit circle. Cauchy's the- 
orem tells that if we know additional integrated quantity (gj\P<j) + f + ) with a set of known functions 
{gi(z),i = l,...,N} one can make the bound stronger as 



det 



X 



\ { gN \p$ + U 



(p$+f+\gl) 

(8n\8i) 



(P$+f + \g N }\ 

(gl\gN) 
(gN\gN) J 



>0. 



(7.8) 



Choosing g n (z) = 7 _ z ^ > Lellouch Q5 8[ ] obtained stronger form factor bounds with statistical anal- 
ysis. Fukunaga and Onogi [ |59| ] improved the bound using also the experimental q 2 spectrum as 
additional inputs. Arnesen et al. [ ]62| ] set g n (z) = z" so that they can obtain the bound on the 
coefficients of the polynomial parameterization of the form factor (z)f(z) = L«=o a nZ n as 



£ kl 2 <x- 



(7.9) 



n=0 



This lead to a great simplification of the problem, although in practice one should truncate the 
polynomial at finite order so that the one has take into account this truncation error as the systematic 
error. Becher and Hill J60|], [ |6lT ] improved Arnesen et al's approach by imposing HQET power 
counting to give stronger constraint than unitarity. Assuming that this power counting argument 
correct they showed that only a few degrees in polynomial is sufficient to approximate the form 
factor. This statement is so far consistent with the observation from the Babar's data in Fig. 10. Of 
course one has to bear in mind that with finite set of data one cannot always exclude the possibility 
that the q 2 spectrum ( z spectrum ) has yet unobserved wiggly behavior from higher order terms in 
the polynomial beyond our experimental resolution, but it will become more clear as experimental 
data will increase. 
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~ E H 2 < 0((A/m h f) Becher-Hill's bound from HQET counting (7.10) 

* n=0 



-e- 



-0.2 





-z 



0.2 



Figure 10: Plot of form factor / + extracted from BaBar experimental data multiplied by a function P<p as a 
function of z- It seems to be consistent with almost linear behavior in z. Figures from [pl|. 



Fermilab collaboration is carrying out an analysis based on Becher-Hill's idea [53]. 

7.2 all-to-all propagators for heavy-light meson 

TrinLat collaborations proposed to construct all-to-all propagators by combining low mode 



averaging [p6[], [|22J and random noise vector technique. The noise should be diluted in time, 
spin and color sources. They have shown that their all-to-all propagator is particularly useful for 
the heavy-light propagator with 20 eigen modes and single random source per dilution for each 
configuration. This method seems very promising. More experience in large volume is needed. 



8. Summary 

Experimental data are offering us a chance to overconstrain CKM. Basic quantities such as 
decay constant, the bag parameter, form factors , b quark masses are important in many ways. 
Several different heavy quark formalism are useful for precision calculation are studied. New 
theoretical or calculation methods are proposed to give better accuracy. 
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